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Abstract — Adaptive OFDMA has recently been recognized as 
a promising technique for providing high spectral efficiency in 
future broadband wireless systems. The research over the last 
decade on adaptive OFDMA systems has focused on adapting 
the allocation of radio resources, such as subcarriers and power, 
to the instantaneous channel conditions of all users. However, 
such "fast" adaptation requires high computational complexity 
and excessive signaling overhead. This hinders the deployment of 
adaptive OFDMA systems worldwide. This paper proposes a slow 
adaptive OFDMA scheme, in which the subcarrier allocation is 
updated on a much slower timescale than that of the fluctuation 
of instantaneous channel conditions. Meanwhile, the data rate 
requirements of individual users are accommodated on the fast 
timescale with high probability, thereby meeting the requirements 
except occasional outage. Such an objective has a natural 
chance constrained programming formulation, which is known 
to be intractable. To circumvent this difficulty, we formulate 
safe tractable constraints for the problem based on recent 
advances in chance constrained programming. We then develop 
a polynomial-time algorithm for computing an optimal solution 
to the reformulated problem. Our results show that the proposed 
slow adaptation scheme drastically reduces both computational 
cost and control signaling overhead when compared with the 
conventional fast adaptive OFDMA. Our work can be viewed as 
an initial attempt to apply the chance constrained programming 
methodology to wireless system designs. Given that most wireless 
systems can tolerate an occasional dip in the quality of service, we 
hope that the proposed methodology will find further applications 
in wireless communications. 

Index Terms — Dynamic Resource Allocation, Adaptive 
OFDMA, Stochastic Programming, Chance Constrained 
Programming 
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I. Introduction 

FUTURE wireless systems will face a growing demand 
for broadband and multimedia services. Orthogonal fre- 
quency division multiplexing (OFDM) is a leading technology 
to meet this demand due to its ability to mitigate wire- 
less channel impairments. The inherent multicarrier nature of 
OFDM facilitates flexible use of subcarriers to significantly 
enhance system capacity. Adaptive subcarrier allocation, re- 
cently referred to as adaptive orthogonal frequency division 
multiple access (OFDMA) 111, Q, has been considered as a 
primary contender in next-generation wireless standards, such 
as 1EEE802.16 WiMAX |3| and 3GPP-LTE |4|. 

In the existing literature, adaptive OFDMA exploits time, 
frequency, and multiuser diversity by quickly adapting sub- 
carrier allocation (SCA) to the instantaneous channel state 
information (CSl) of all users. Such "fast" adaptation suffers 
from high computational complexity, since an optimization 
problem required for adaptation has to be solved by the base 
station (BS) every time the channel changes. Considering the 
fact that wireless channel fading can vary quickly (e.g., at 
the order of milli-seconds in wireless cellular system), the 
implementation of fast adaptive OFDMA becomes infeasible 
for practical systems, even when the number of users is small. 
Recent work on reducing complexity of fast adaptive OFDMA 
includes 15], etc. Moreover, fast adaptive OFDMA requires 
frequent signaling between the BS and mobile users in order 
to inform the users of their latest allocation decisions. The 
overhead thus incurred is likely to negate the performance 
gain obtained by the fast adaptation schemes. To date, high 
computational cost and high control signaling overhead are 
the major hurdles that prevent adaptive OFDMA from being 
deployed in practical systems. 

We consider a slow adaptive OFDMA scheme, which is 
motivated by [7|, to address the aforementioned problem. 
In contrast to the common belief that radio resource allo- 
cation should be readapted once the instantaneous channel 
conditions change, the proposed scheme updates the SCA 
on a much slower timescale than that of channel fluctuation. 
Specifically, the allocation decisions are fixed for the duration 
of an adaptation window, which spans the length of many 
coherence times. By doing so, computational cost and control 
signaling overhead can be dramatically reduced. However, this 
implies that channel conditions over the adaptation window are 
uncertain at the decision time, thus presenting a new challenge 
in the design of slow adaptive OFDMA schemes. An important 
question is how to find a valid allocation decision that remains 
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optimal and feasible for the entire adaptation window. Such 
a problem can be formulated as a stochastic programming 
problem, where the channel coefficients are random rather than 
deterministic. 

Slow adaptation schemes have recently been studied in other 
contexts such as slow rate adaptation [7J, [8J and slow power 
allocation Therein, adaptation decisions are made solely 
based on the long-term average channel conditions instead of 
fast channel fading. Specifically, random channel parameters 
are replaced by their mean values, resulting in a deterministic 
rather than stochastic optimization problem. By doing so, 
quality-of-service (QoS) can only be guaranteed in a long-term 
average sense, since the short-term fluctuation of the channel is 
not considered in the problem formulation. With the increasing 
popularity of wireless multimedia applications, however, there 
will be more and more inelastic traffic that require a guarantee 
on the minimum short-term data rate. As such, slow adaptation 
schemes based on average channel conditions cannot provide 
a satisfactory QoS. 

On another front, robust optimization methodology can be 
applied to meet the short-term QoS. For example, robust opti- 
mization method was applied in llQl- lfTTl to find a solution that 
is feasible for the entire uncertainty set of channel conditions, 
i.e., to guarantee the instantaneous data rate requirements 
regardless of the channel realization. Needless to say, the 
resource allocation solutions obtained via such an approach 
are overly conservative. In practice, the worst-case channel 
gain can approach zero in deep fading, and thus the resource 
allocation problem can easily become infeasible. Even if the 
problem is feasible, the resource utilization is inefficient as 
most system resources must be dedicated to provide guarantees 
for the worst-case scenarios. 

Fortunately, most inelastic traffic such as that from mul- 
timedia applications can tolerate an occasional dip in the 
instantaneous data rate without compromising QoS. This 
presents an opportunity to enhance the system performance. 
In particular, we employ chance constrained programming 
techniques by imposing probabilistic constraints on user QoS. 
Although this formulation captures the essence of the problem, 
chance constrained programs are known to be computationally 
intractable except for a few special cases fT2l|. In general, such 
programs are difficult to solve as their feasible sets are often 
non-convex. In fact, finding feasible solutions to a generic 
chance constrained program is itself a challenging research 
problem in the Operations Research community. It is partly 
due to this reason that the chance constrained programming 
methodology is seldom pursued in the design of wireless 
systems. 

In this paper, we propose a slow adaptive OFDMA scheme 
that aims at maximizing the long-term system throughput 
while satisfying with high probability the short-term data 
rate requirements. The key contributions of this paper are as 
follows: 

• We design the slow adaptive OFDMA system based on 
chance constrained programming techniques. Our formu- 
lation guarantees the short-term data rate requirements 
of individual users except in rare occasions. To the best 
of our knowledge, this is the first work that uses chance 



constrained programming in the context of resource allo- 
cation in wireless systems. 

• We exploit the special structure of the probabilistic 
constraints in our problem to construct safe tractable 
constraints (STC) based on recent advances in the chance 
constrained programming literature. 

• We design an interior-point algorithm that is tailored for 
the slow adaptive OFDMA problem, since the formu- 
lation with STC, although convex, cannot be trivially 
solved using off-the-shelf optimization software. Our 
algorithm can efficiently compute an optimal solution to 
the problem with STC in polynomial time. 

The rest of the paper is organized as follows. In Section 
nil we discuss the system model and problem formulation. An 
STC is introduced in Section [III] to solve the original chance 
constrained program. An efficient tailor-made algorithm for 
solving the approximate problem is then proposed in Section 
IrVl In Section [V] we reduce the problem size based on some 
practical assumptions, and show that the revised problem can 
be solved by the proposed algorithm with much lower com- 
plexity. In Section [Vll the performance of the slow adaptive 
OFDMA system is investigated through extensive simulations. 
Finally, the paper is concluded in Section IVIII 



II. System Model and Problem Formulation 

This paper considers a single-cell multiuser OFDM system 
with K users and N subcarriers. We assume that the instan- 
taneous channel coefficients of user k and subcarrier n are 
described by complex GaussiarQ random variables h^^\ ~ 
CA/'(0,CT^), independent in both 11 and k. The parameter Uk 
can be used to model the long-term average channel gain as 

<yk — \ ^] 'Sk, where dk is the distance between the BS and 



subscriber k, do is the reference distance, 7 is the amplitude 
path-loss exponent and Sk characterizes the shadowing effect. 



Hence, the channel gain g[*^^ — l^k'nl" ™ exponential 
random variable with probability density function (PDF) given 
by 



t) 



AO 



■ exp 



(1) 



The transmission rate of user k on subcarrier n at time t is 
given by 



Wlog^ 1 



(t) ^ 

Pt9kjn 
TNn 



where pt is the transmission power of a subcarrier, g^*^^^ is the 
channel gain at time t, W is the bandwidth of a subcarrier, A^o 
is the power spectral density of Gaussian noise, and F is the 
capacity gap that is related to the target bit error rate (BER) 
and coding-modulation schemes. 

In traditional fast adaptive OFDMA systems, SCA decisions 
are made based on instantaneous channel conditions in order 

'Although the techniques used in this paper are applicable to any fading 
distribution, we shall prescribe to a particular distribution of fading channels 
for illustrative purposes. 

-The case when frequency con'elations exist among subcarriers will be 
discussed in Section IvTl 
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Fig. 1. Adaptation timescales of fast and slow adaptive OFDMA system 
(SCA = SubCarrier Allocation). 



to maximize the system throughput. As depicted in Fig. [T^, 
SCA is performed at the beginning of each time slot, where 
the duration of the slot is no larger than the coherence time of 
the channel. Denoting by x"^\ the fraction of airtime assigned 
to user k on subcarrier n, fast adaptive OFDMA solves at each 
time slot t the following linear programming problem: 



max 

„(') 



s.t. 



K N 

EE 



At) Jt) 

^k,n' k,n 



k=l n=l 

N 

E(*) (t) ^ 

Tl=l 

E4!l<i, 

k=l 

(t) 
k.n 



(2) 



(3) 



X 



Vn 



> 0, Vfc, n, 



where the objective function in (|2]l represents the total system 
throughput at time t, and (O represents the data rate constraint 
of user k at time t with qk denoting the minimum required 
data rate. We assume that is known by the BS and can be 
different for each user k. Since g^l^ (and hence r^*^^) varies on 
the order of coherence time, one has to solve the Problem Vhst 
at the beginning of every time slot t to obtain SCA decisions. 
Thus, the above fast adaptive OFDMA scheme is extremely 
costly in practice. 

In contrast to fast adaptation schemes, we propose a slow 
adaptation scheme in which SCA is updated only every adap- 
tation window of length T. More precisely, SCA decision is 
made at the beginning of each adaptation window as depicted 
in Fig. [TJ), and the allocation remains unchanged till the next 
window. We consider the duration T of a window to be 
large compared with that of fast fading fluctuation so that 
the channel fading process over the window is ergodic; but 
small compared with the large-scale channel variation so that 
path-loss and shadowing are considered to be fixed in each 
window. Unlike fast adaptive systems that require the exact 
CSI to perform SCA, slow adaptive OFDMA systems rely 
only on the distributional information of channel fading and 
make an SCA decision for each window. 

Let Xk,n € [0, 1] denote the SCA for a given adaptation 



windo\\0. Then, the time-average throughput of user k during 
the window becomes 



N 

^k ^ ^ '^k^n'^kjU-; 
n=l 



where 



1 



J .It 



is the time-average data rate of user k on subcarrier n during 
the adaptation window. The time-average system throughput 
is given by 

K N 
fc=l n=l 

Now, suppose that each user has a short-term data rate require- 
ment qk defined on each time slot. If X]^=i ^k,n'r^k \ < ^k, 
then we say that a rate outage occurs for user k at time slot t, 
and the probability of rate outage for user k during the window 
[io, + T] is defined as 



P^°"' ^ Pr 



' N 

Xk^nrf^^ < qk 



Vt G [io,io + T], 



where to is the beginning time of the window. 

Inelastic applications, such as voice and multimedia, that 
are concerned with short-term QoS can often tolerate an 
occasional dip in the instantaneous data rate. In fact, most 
applications can run smoothly as long as the short-term data 
rate requirement is satisfied with sufficiently high probability. 
With the above considerations, we formulate the slow adaptive 
OFDMA problem as follows: 

K N 

T'siow: max E ^^''^"'^ {'^^n} 



k=l n=l 
( N 



.t. Pr ^ 



K 



Xk 



n' k.n 



>qk>>l-ek, Vfc (5) 



,Xk,n < 1, Vn 
k=l 

Xk.n > 0, Vfc,n, 

where the expectatior0 in (|4|i is taken over the random channel 
process g = {g^*^} for t G [to, to + T], and efc G [0, 1] in © 
is the maximum outage probability user k can tolerate. In the 
above formulation, we seek the optimal SCA that maximizes 
the expected system throughout while satisfying each user's 
short-term QoS requirement, i.e., the instantaneous data rate 
of user k is higher than qk with probability at least 1 — e^. 
The above formulation is a chance constrained program since 
a probabilistic constraint (|5]) has been imposed. 



^^It is practical to assume „ as a real number in slow adaptive OFDMA. 
Since the data transmitted during each window consists of a large mount of 
OFDM symbols, the time-sharing factor „ can be mapped into the ratio 
of OFDM symbols assigned to user k for transmission on subcarrier n. 

^In (4), we replace the time-average data rate fj. „ by its ensemble average 
E l*"^*^! due to the ergodicity of channel fading over the window. 
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III. Safe Tractable Constraints 



Despite its utility and relevance to real applications, the 
chance constraint ^ imposed in T'siow makes the optimization 
highly intractable. The main reason is that the convexity of the 
feasible set defined by ^ is difficult to verify. Indeed, given 
a generic chance constraint Pr {F{x, r) > 0} < e where r is a 
random vector, x is the vector of decision variable, and F is a 
real- valued function, its feasible set is often non-convex except 

for very few special cases lfT2l . ifTSl . Moreover, even with the 

(t) 



Then, the allocation decision x satisfies 

Pr |^^Xfc,„ri*^ > gfcj > 1 - efc, Vfc. (12) 

Proof: Our argument will use the Bernstein approxi- 
mation theorem proposed in Suppose there exists an 
X e M^^'' such that Gfe(x) < 0, i.e.. 



N 



inf l^fe + Afe(-£i ^ifc,n) - flogefcj < 0. (13) 



qk + logE <^ exp f - 

n=l ^ 



nice function in (|5]l, i.e., F{x,r) = qu — Sn=i ^fe.n^i i The function inside the infg>o{-} is equal to 
bilinear in x and r, with independent entries r^'^ in r whose 
distribution is known, it is still unclear how to compute the 
probability in (|5]l efficiently. 

To circumvent the above hurdles, we propose the following 
formulation T^siow by replacing the chance constraints (|5]) with 
a system of constraints Ti. such that (i) x is feasible for (|5]l 
whenever it is feasible for T-L, and (ii) the constraints in T-L are 
convex and efficiently computablqfj. The new formulation is 
given as follows: 



^ik..4\]\ - Q\ogek (14) 



= qk + g log E { cxp Q 



n=l ^ J 



(15) 



N 



£>logE ^ exp( g ^qk-^l ^fc."*"*," ) )} ^ fi'logefc, 




K N 



Vsiow ■ max Xk^n^ {'"£1} 

k=ln=l 

s.t. mi<qk + gy^Aki-g ^Xk,n) 

n=l 



(16) 



(6) 



K 

^Xk,n < 1, Vn 

k=l 

Xk.n > 0, V/c, n, 



glogek}<0, yk (7) 

(8) 
(9) 



where A/c(-) is the cumulant generating function of r[*^ 



Afc(-e ife,„) = log 



1 



1 

• — exp 



(10) 



In the following, we first prove that any solution x that is 
feasible for the STC (|7]i in "Psiow is also feasible for the chance 
constraints Q. Then, we prove that T^siow is convex. 

Proposition 1. Suppose that g[.*^ (and hence r^k\j 
are independent random variables for different n and 
k, where the PDF of g^^^ follows ([T]l. Furthermore, 
given Efc > 0, suppose that there exists an ii. = 



[a;i_i,--- ,xn,1t ■ ■ ,xi^Kt - ■ ■ ,xn,k\ 



pNK 



such that 



Gfc(x) = inf <gA; + £)y^ Afe(-g ^Xfc,.„)-£»logefe> <0, Vfc. 

(11) 



where the expectation E {•} can be computed using the distri- 
butional information of g^'^ in ([TJ, and ( fT5] l follows from the 
independence of random variable r^*|j^ over n. 

Let J^fc(x, r) = Qk - J2n=i ^k,nr'j^^n- Then, (O is equiva- 
lent to 

M{gE{cxp{g-^Fk{S^,r))}^gek} < 0. (17) 

According to Theorem |2] in Appendix |A] the chance con- 
straints (fT2] l hold if there exists a g > satisfying (fTTI i. Thus, 
the validity of ( fTSl i is guaranteed by the validity of (fTTT i. □ 
Now, we prove the convexity of (|7]i in the following 
proposition. 

Proposition 2. The constraints imposed in ^ are convex in 
X = [xia, • • • , XN,i, • ■ ■ , xi^K, • • • , xn,k]'^ e R^^. 

Proof: The convexity of (Q can be verified as follows. 
We define the function inside the infg>o{-} in ( fTTI ) as 

N 

Hkix, g) ^ qk + g^Aki-g^^Xk,n) - g^ogek, Vfc. (18) 

n=l 

It is easy to verify the convexity of Hk{x, g) in (x, g), since 
the cumulant generating function is convex. Hence, Gfe(x) in 
( fTTI ) is convex in x due to the preservation of convexity by 
minimization over g > 0. □ 



IV. Algorithm 

In this section, we propose an algorithm for solving Problem 
T^siow In Alow, the STC O arises as a subproblem, which 
by itself requires a minimization over g. Hence, despite its 
convexity, the entire problem Psiow cannot be trivially solved 



^Condition (i) is refen'ed to as "safe" condition, and condition (ii) is referred 
to as "tractable" condition. 



^For the reader's convenience, botli tlie theorem and a rough proof are 
provided in Appendix IaI 
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Algorithm 1 Structure of the Proposed Algorithm 

Require: The feasible solution set of Problem T'siow is a 

compact set X defined by (l7]i-(|9l). 
1: Construct a polytope X° D A" by ((HI)-©. Set i 0. 
2: Choose a query point {Subsection |73A-i) at the ith 

iteration as x' by computing the analytic center of X'. 

Initially, set x'' = e/if G where e is an iV- vector of 

ones. 

3: Query the separation oracle (Subsection \rV\A-I) with x*: 
4: if X* e A" then 

5: generate a hyperplane (optimality cut) through x* to 
remove the part of X"^ that has lower objective values 
6: else 

7: generate a hyperplane (feasibility cut) through x* to 
remove the part of X"^ that contains infeasible solutions. 
8: end if 

9: Set i -i— i + 1, and update AT'^^ by the separation 
hyperplane. 

10: if termination criterion {Subsection^^B) is satisfied then 
II: stop 
12: else 

13: return to step 2. 
14: end if 



using standard solvers of convex optimization. This is due 
to the fact that the subproblem introduces difficulties, for 
example, in defining the barrier function in path-following 
algorithms or providing the (sub-)gradient in primal-dual 
methods (see flA'l for details of these algorithms). Fortunately, 
we can employ interior point cutting plane methods to solve 
Problem T'siow (see lITSi for a survey). Before we delve into 
the details, let us briefly sketch the principles of the algorithm 
as follows. 

Suppose that we would like to find a point x that is feasible 
for ©-(IHl and is within a distance of i5 > to an optimal 
solution X* of "Psiow, where (5 > is an error tolerance 
parameter (i.e., x satisfies ||x — x*||2 < 5). We maintain the 
invariant that at the beginning of each iteration, the feasible 
set is contained in some polytope (i.e., a bounded polyhedron). 
Then, we generate a query point inside the polytope and ask 
a "separation oracle" whether the query point belongs to the 
feasible set. If not, then the separation oracle will generate 
a so-called separating hyperplane through the query point to 
cut out the polytope, so that the remaining polytope contains 
the feasible setlj Otherwise, the separation oracle will return 
a hyperplane through the query point to cut out the polytope 
towards the opposite direction of improving objective values. 

We can then proceed to the next iteration with the new 
polytope. To keep track of the progress, we can use the so- 
called potential value of the polytope. Roughly speaking, when 
the potential value becomes large, the polytope containing the 
feasible set has become small. Thus, if the potential value 
exceeds a certain threshold, so that the polytope is negligibly 
small, then we can terminate the algorithm. As will be shown 



later, such an algorithm will in fact terminate in a polynomial 
number of steps. 

We now give the structure of the algorithm. A detailed flow 
chart is shown in Fig. |2]for readers' interest. 



Initialize: 
The set X'^ : {A° ,h°) 
and xO = e/K. 



-> e*=arg inf [H(x',£i)] 

Q>0 



Oracle 




Adding 
feasibility cut 
(23) 




Adding 
optimality cut 
(24) 




Update 
' X^+l:(A^+l,b^+l) 




* 


* 
















The problem is 
infeasible. 




The problem is feasible. 
The optimal solution x* = x' . 











End 



Fig. 2. Flow chart of the algorithm for solving Problem Pslow 



A. The Cutting-Plane-Based Algorithm 

1 ) Query Point Generator: ( Step 2 in Algorithm Q} 
In each iteration, we need to generate a query point inside 
the polytope A*. For algorithmic efficiency, we adopt the ana- 
lytic center (AC) of the containing polytope as the query point 
|17|. The AC of the polytope A* = {x e R^^ : yl*x < b*} 
at the ith iteration is the unique solution x* to the following 
convex problem: 



'Note that such a separating hyperplane exists due to the convexity of the 
feasible set 1161 . 



M' 

max log sl^ 

s.t. = b* - 



(19) 
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We define the optimal value of the above problem as the 
potential value of the poly tope X"^. Note that the uniqueness 
of the analytic center is guaranteed by the strong convexity 
of the potential function s* — X)m=i ^^S ■^m' assuming that 
X* is bounded and has a non-empty interior The AC of a 
polytope can be viewed as an approximation to the geometric 
center of the polytope, and thus any hyperplane through the 
AC will separate the polytope into two parts with roughly the 
same volume. 

Although it is computationally involved to directly solve 
( fT9] l in each iteration, it is shown in fT8 1 that an approximate 
AC is sufficient for our purposes, and that an approximate AC 
for the (i + l)st iteration can be obtained from an approximate 
AC for the ith iteration by applying 0(1) Newton steps. 
2) The Separation Oracle: (Steps 3-8 in Algorithm\I} 
The oracle is a major component of the algorithm that plays 
two roles: checking the feasibility of the query point, and 
generating cutting planes to cut the current set. 

• Feasibility Check 

We write the constraints of Psiow in a condensed form as 
follows: 



Gfe(x) = inf {Hki^,g)}<0, 
A° X < b" 



Vfc 



(20) 
(21) 



where 



A 



In In • ■ ■ In 

—Ink 



,,(N+NK)xNK 



pN+NK 



with In and gat denoting the N x N identity matrix and A^- 
vector of ones respectively, and (l2Tl i is the combinatior0 of ([8]l 
and (|9]l. Now, we first use (l2Tl i to construct a relaxed feasible 
set via 

X° = {x e R^^ : A° X < b"}. (22) 

Given a query point x e X", we can verify its feasibility to 
^siow by checking if it satisfies (|20] |. i.e., if infg>o{-fffe(x, g)} 
is no larger than 0. This requires solving a minimization 
problem over g > 0. Due to the unimodality of g) in g, 

we can simply take a line search procedure, e.g., using Golden- 
section search or Fibonacci search, to find the minimizer 
g*. The line search is more efficient when compared with 
derivative-based algorithms, since only function evaluation^ 
are needed during the search. 

• Cutting Plane Generation 

In each iteration, we generate a cutting plane, i.e., a hyper- 
plane through the query point, and add it as an additional 
constraint to the current polytope X*. By adding cutting 
plane(s) in each iteration, the size of the polytope keeps shrink- 
ing. There are two types of cutting planes in the algorithm 
depending on the feasibility of the query point. 

If the query point x* £ is infeasible, then a hyperplane 

*To reduce num erical errors in computation, we suggest normalizing each 
constraint in j2H . 

'The cumulant generating function Afe(-) in jlOt can be evaluated numer- 
ically, e.g., using rectangular rule, trapezoid rule, or Simpson's rule, etc. 



csMqA feasibility cut is generated at x* as follows: 

T 

(x-x')<o, yaeK, 



(23) 



where || • || is the Euclidean norm, K = {k : 
iJfc(x%r) > 0, fc = 1,2,--- ,7^} is the set of 
users whose chance constraints are violated, and u' " = 



' ' ' ' ^Ar,i' ■ • ■ ' ' ' ' i '^n.kI 
ent of Gk(x) with respect to x, i.e., 

t,R _ <9i^s(x,^?*) 



fNK 



is the gradi- 



dxi 



w 

In 2 Jo 



Pit. 

TNo 



In 1 



Mi. 

FNo I a. 



■ exp 



The reason we call ( l23T l a feasibility cut(s) is that any x which 
does not satisfy (l23b must be infeasible and can hence be 
dropped. 

If the point x' is feasible, then an optimality cut is generated 
as follows: „ 

V 



where v= [ - E{r['l} , ■ ■ ■ , -^{4],} , 



(x - X*) < 0, 



(24) 



fNK 



is the derivative of the objective of 



T-'.slow in I© with respect to x. The reason we call (l24l i an 
optimality cut is that any optimal solution x* must satisfy (l24l) 
and hence any x which does not satisfy (l24l) can be dropped. 

Once a cutting plane is generated according to ( |23] ) or ( l24b . 
we use it to update the polytope X* at the ith iteration as 
follows 

X' = {xe M^-^' : A'x < W}. 

Here, A* and b* are obtained by adding the cutting plane to the 
previous polytope X*^^. Specifically, if the oracle provides a 
feasibility cut as in ( |23] ), then 



A' = 
b' = 



1,1-1 



u. 



^(M'-^ + \K\)xNK 



where Mi_i is the number of rows in Ai^i, and | • | is the 
number of elements contained in the given set; if the oracle 
provides an optimality cut as in (|24|) . then 



A' 



(v/||v|D 
(v 



p(A/*"Hl)xAfif 



vin^x' 



B. Global Convergence & Complexity (Step 10 in AlgorithmU^ 

In the following, we investigate the convergence properties 
of the proposed algorithm. As mentioned earlier, when the 
polytope is too small to contain a full-dimensional closed ball 
of radius 6 > 0, the potential value will exceeds a certain 
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threshold. Then, the algorithm can terminate since the query 
point is within a distance of (5 > to some optimal solution of 
Vsiow Such an idea is formalized in |18 |, where it was shown 
that the analytic center-based cutting plane method can be used 
to solve convex programming problems in polynomial time. 
Upon following the proof in l,18J . we obtain the following 
result: 

Theorem 1. (cf. ^7g|/ ) Let 5 > be the error tolerance 
parameter, and let m be the number of variables. Then, 
Algorithm [7] terminates with a solution x that is feasible for 
Vsiow ond satisfies ||x — x*||2 < S for some optimal solution 
X* to Vsiow after at most 0{{m/ 5)"^) iterations. 

Thus, the proposed algorithm can solve Problem T^siow 
within 0{{N K / SY) iterations. It turns out that the algo- 
rithm can be made considerably more efficient by dropping 
constraints that are deemed "unimportant" in |T9l. By incor- 
porating such a strategy in Algorithm [T] the total number 
of iterations needed by the algorithm can be reduced to 
0{NK\og^{l/5)). We refer the readers to US, Ull for 
details. 

C. Complexity Comparison between Slow and Fast Adaptive 
OFDMA 

It is interesting to compare the complexity of slow and 
fast adaptive OFDMA schemes formulated in T'siow and T'fast, 
respectively. To obtain an optimal solution to T'fast, we need 
to solve a linear program (LP). This requires ©(V NKLq) 
iterations, where Lq is number of bits to store the data defining 
the LP II20I . At first glance, the iteration complexity of solving 
a fast adaptation Vhst can be lower than that of solving T'siow 
when the number of users or subcarriers are large. However, it 
should be noted that only one "Psiow needs to be solved for each 
adaptation window, while Vf^st has to be solved for each time 
slot. Since the length of adaptation window is equal to T time 
slots, the overall complexity of the slow adaptive OFDMA 
can be much lower than that of conventional fast adaptation 
schemes, especially when T is large. 

Before leaving this section, we emphasize that the advantage 
of slow adaptive OFDMA lies not only in computational cost 
reduction, but also in reducing control signaling overhead. We 
will investigate this in more detail in Section IVll 

V. Problem Size Reduction 

In this section, we show that the problem size of T^siow can be 
reduced from NK variables to K variables under some mild 
assumptions. Consequently, the computational complexity of 
slow adaptive OFDMA can be markedly lower than that of 
fast adaptive OFDMA. 

In practical multicarrier systems, the frequency intervals 
between any two subcarriers are much smaller than the carrier 
frequency. The reflection, refraction and diffusion of electro- 
magnetic waves behave the same across the subcarriers. This 
implies that the channel gain g'j^]^ is identically distributed 
over n (subcarriers), although it is not needed in our algorithm 
derivations in the previous sections. 



When g^*^ for different n are identically distributed, dif- 
ferent subcarriers become indistinguishable to a user k. In 
this case, the optimal solution, if exists, does not depend on 
n. Replacing Xk.n by Xk in T^siow, we obtain the following 
formulation: 

K N 

'^siow^ max 

^ fc=ln=l 

s.t. mi{qk + gNAk{-g^^Xk)~g\ogek}<0, Vk 

g>0 
K 

fc=l 

Xk > 0, Vfc. 

Note that the problem structure of V^^,^^ is exactly the 
same as that of T-'siow, except that the problem size is re- 
duced from NK variables to K variables. Hence, the al- 
gorithm developed in Section |IV] can also be applied to 
solve "PjIo^j,, with the following vector/matrix size reductions: 
= [bn^-IkV e M(i+^)x^, bO = [1, o,--- , 0]^ e 
in (ED, u*'"" = WT.--- ,uj^V e in (|23]l, and 
V = [ - E{rf^}, • • • , -E{r^^}]^ € in (El. Compared 
with Psiow, the iteration complexity of P^io^ is now reduced 
to 0{K\og^{l/5)). Indeed, this can even be lower than the 
complexity of solving one T'fast — 0{\/ NKLq), since K 
is typically much smaller than N in real systems. Thus, the 
overall complexity of slow adaptive OFDMA is significantly 
lower than that of fast adaptation over T time slots. 

Before leaving this section, we emphasize that the problem 
size reduction in Vl^^^ does not compromise the optimality of 
the solution. On the other hand, T'siow is more general in the 
sense that it can be applied to systems in which the frequency 
bands of parallel subchannels are far apart, so that the channel 
distributions are not identical across different subchannels. 

VI. Simulation Results 

In this section, we demonstrate the performance of our 
proposed slow adaptive OFDMA scheme through numerical 
simulations. We simulate an OFDMA system with 4 users 
and 64 subcarriers. Each user k has a requirement on its 
short-term data rate qk = 20bps. The 4 users are assumed 
to be uniformly distributed in a cell of radius R — 100m. 
That is, the distance dk between user k and the BS follows 
the distributioi{3 f{d) = j^. The path-loss exponent 7 is 
equal to 4, and the shadowing effect Sk follows a log-normal 
distribution, i.e., 10 log]^o('Sfe) ^ A/'(0, 8dB). The small-scale 
channel fading is assumed to be Rayleigh distributed. Suppose 
that the transmission power of the BS on each subcarrier is 
90dB measured at a reference point 1 meter away from the 
BS, which leads to an average received power of lOdB at the 
boundary of the celj*^ In addition, we set W = IHz and A^o = 
1, and the capacity gap is L = - log(5BER)/1.5 = 5.0673, 

'"The distribution of user's distance from tlie BS f{d) = is derived 
from the uniform distribution of user's position f(x,y) = ^j^^ , where (x, y) 
is the Cartesian coordinate of the position. 

"The average received power at the boundary is calculated by 90dB + 
lOlogio {^)^^dB = lOdB due to the path-loss effecL 
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5 10 15 20 25 30 10 20 30 40 50 60 

Iteration Window Number 



Fig. 3. Trace of the difference of objective value fe' between adjacent 
iterations {e^ = 0.2). 



Fig. 4. Number of iterations for convergence of all the feasible windows 
= 0.2). 



where the target BER is set to be 10^^. Moreover, the length 
of one slot, within which the channel gain remains unchanged, 
IS To = lms@ The length of the adaptation window is chosen 
to be T = Is, implying that each window contains 1000 
slots. Suppose that the path loss and shadowing do not change 
within a window, but varies independently from one window to 
another For each window, we solve the size-reduced problem 
V!.fg^, and later Monte-Carlo simulation is conducted over 61 
independent windows that yield non-empty feasible sets of 
n'low when efe = 0.1. 

In Fig. [3] and Fig. |4] we investigate the fast convergence 
of the proposed algorithm. The error tolerance parameter is 
chosen as 5 — 10^^. In Fig. |3] we record the trace of one 
adaptation windowT^ and plot the improvement in the objective 
function value (i.e., system throughput) in each iteration, i.e., 
Ab — If — . When A& is positive, the objective value 
increases with each iteration. It can be seen that A6 quickly 
converges to close to zero within only 27 iterations. We also 
notice that fluctuation exists in A& within the first 1 1 iterations. 
This is mainly because during the search for an optimal 
solution, it is possible for query points to become infeasible. 
However, the feasibility cuts (|23] | then adopted will make sure 
that the query points in subsequent iterations will eventually 
become feasible. The curve in Fig. [3] verifies the tendency. 
As T^siow is convex, this observation implies that the proposed 
algorithm can converge to an optimal solution of "Psiow within 
a small number of iterations. In Fig. 2] we plot the number 
of iterations needed for convergence for different appUcation 
windows. The result shows that the proposed algorithm can 
in general converge to an optimal solution of "Psiow within 
35 iterations. On average, the algorithm converges after 22 



12' 



The coherence time is given by To 



9e 



where c is the speed of 



167r/c V ' 

light, fc is the carrier frequency, and v is the velocity of mobile user As an 
example, we choose fc = 2.5GHz, and if the user is moving at 45 miles per 
hour, the coherence time is around 1ms. 

"The simulation results show that all the feasible windows appear with 
similar convergence behavior 



X feasibie 
o infeasibie 




10 20 30 40 50 60 70 80 90 100 
Window Number 



Fig. 5. Number of iterations for feasibility check of all the windows (ej. : 
0.2). 



iterations, where each iteration takes 1.467 seconds!*^ 

Moreover, we plot the number of iterations needed for 
checking the feasibility of T'siow In Fig. |5] we conduct a 
simulation over 100 windows, which consists of 61 feasible 
windows (dots with cross) and 39 infeasible windows (dots 
with circle). On average, the algorithm can determine if "Psiow 
is feasible or not after 7 iterations. The quick feasibility check 
can help to deal with the admission of mobile users in the 
cell. Particularly, if there is a new user moving into the cell, 
the BS can adopt the feasibility check to quickly determine 
if the radio resources can accommodate the new user without 
sacrificing the current users' QoS requirements. 

In Fig. |6] we compare the spectral efficiency of slow adap- 

'''We conduct a simulation on Matlab 7.0.1, where the system configu- 
rations are given as: Processor: Intel(R) Core(TM)2 CPU P8400@2.26GHz 
2.27GHz, Memory: 2.00GB, System Type: 32-bit Operating System. 
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- fast adaptation 

- slow adaptation (e^. = 0. 1) 




outage probability of user 1 



20 30 40 

Window Number 




0.1 

0.05 



20 30 40 

outage probability of user 4 



£i. = 0.3 



Fig. 6. Comparison of system spectral efficiency between fast adaptive Fig. 7. Outage probability of the 4 users over 61 independent feasible 
OFDMA and slow adaptive OFDMA. windows. 



tive OFDMA with that of fast adaptive OFDMA^l, where zero 
outage of short-term data rate requirement is ensured for each 
user. In addition, we take into account the control overheads 
for subcarrier allocation, which will considerably affect the 
system throughput as well. Here, we assume that the control 
signaling overhead consumes a bandwidth equivalent to 10% 
of a slot length Tq every time SCA is updated 1211 . Note that 
within each window that contains 1000 slots, the control sig- 
naling has to be transmitted 1000 times in the fast adaptation 
scheme, but once in the slow adaptation scheme. In Fig. |6] 
the line with circles represents the performance of the fast 
adaptive OFDMA scheme, while that with dots corresponds 
to the slow adaptive OFDMA. The figure shows that although 
slow adaptive OFDMA updates subcarrier allocation 1000 
times less frequently than fast adaptive OFDMA, it can achieve 
on average 71.88% of the spectral efficiency. Considering the 
substantially lower computational complexity and signaling 
overhead, slow adaptive OFDMA holds significant promise 
for deployment in real-world systems. 

As mentioned earlier, VsXovi is more conservative than the 
original problem T^siow, implying that the outage probability is 
guaranteed to be satisfied if subcarriers are allocated according 
to the optimal solution of "Psiow This is illustrated in Fig. |7] 
which shows that the outage probability is always lower than 
the desired threshold = 0.1. 

Fig. |7] shows that the subcarrier allocation via T'siow could 
still be quite conservative, as the actual outage probability is 
much lower than e^. One way to tackle the problem is to set 
€k to be larger than the actual desired value. For example, we 
could tune €k from 0.1 to 0.3. By doing so, one can potentially 
increase the system spectral efficiency, as the feasible set of 

'^For illustrative purpose, we have only considered Vf-^a as one of the 
typical formulations of fast adaptive OFDMA in our comparisons. However, 
we should point out that there are some work on fast adaptive OFDMA which 
impose less restrictive constraints on user data rate requirement. For example, 
in [5J, it considered average user data rate constraints which exploits time 
diversity to achieve higher spectral efficiency. 




Fig. 8. Spectral efficiency versus tolerance parameter e^. Calculated from 
the average overall system throughput on one window, where the long-term 
average channel gain o-fe of the 4 users are — 65.11dB, — 56.28dB, — 68.14dB 
and — 81.96dB, respectively. 



T'slow is enlarged. A question that immediately arises is how 
to choose the right e^, so that the actual outage probability 
stays right below the desired value. Towards that end, we can 
perform a binary search on to find the best parameter that 
satisfies the requirement. Such a search, however, inevitably 
involves high computational costs. On the other hand. Fig. |8] 
shows that the gain in spectral efficiency by increasing is 
marginal. The gain is as little as 0.5 bps/Hz/subcarrier when 
tk is increased drastically from 0.05 to 0.7. Hence, in practice, 
we can simply set to the desired outage probability value 
to guarantee the QoS requirement of users. 

In the development of the STC (jT), we considered that the 
channel gain g^. „ are independent for different n's and fc's. 
While it is true that channel fading is independent across 
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optimization problem. Based on this formulation, we design a 
polynomial-time algorithm for subcarrier allocation in slow 
adaptive OFDMA. Our simulation results showed that the 
proposed algorithm converges within 22 iterations on average. 

In the future, it would be interesting to investigate the 
chance constrained subcarrier allocation problem when fre- 
quency correlation exists, or when the channel distribution 
information is not perfectly known at the BS. Moreover, it is 
worthy to study the tightness of the Bernstein approximation. 
Another interesting direction is to consider discrete data rate 
and exclusive subcarrier allocation. In fact, the proposed 
algorithm based on cutting plane methods can be extended 
to incorporate integer constraints on the variables (see e.g., 

mi 

Finally, our work is an initial attempt to apply the chance 
constrained programming methodology to wireless system 
designs. As probabilistic constraints arise quite naturally in 
many wireless communication systems due to the randomness 
in channel conditions, user locations, etc., we expect that 
chance constrained programming will find further applications 
in the design of high performance wireless systems. 



Fig. 9. Compaiison of outage probability of 4 users with and without 
frequency correlations in channel model. 



Appendix A 
Bernstein Approximation Theorem 



different users, it is typically correlated in the frequency 
domain. We investigate the effect of channel correlation in 
frequency domain through simulations. A wireless channel 
with an exponential decaying power profile is adopted, where 
the root-mean-square delay is equal to 37.79ns. For com- 
parison, the curves of outage probability with and without 
frequency correlation are both plotted in Fig. |9l We choose 
the tolerance parameter to be efe — 0.3. The figure shows 
that with frequency-domain correlation, the outage probability 
requirement of 0.3 is violated occasionally. Intuitively, such 
a problem becomes negligible when the channel is highly 
frequency selective, and is more severe when the channel is 
more frequency flat. To address the problem, we can set ek 
to be lower than the desired outage probability valu43. For 
example, when we choose = 0.1 in Fig. |9] the outage 
probabilities all decreased to lower than the desired value 0.3, 
and hence the QoS requirement is satisfied (see the line with 
dots). 

VII. Conclusions 

This paper proposed a slow adaptive OFDMA scheme 
that can achieve a throughput close to that of fast adaptive 
OFDMA schemes, while significantly reducing the computa- 
tional complexity and control signaling overhead. Our scheme 
can satisfy user data rate requirement with high probability. 
This is achieved by formulating our problem as a stochastic 

Alternatively, we can divide N subcarriers into subchannels (each 
subchannel consists Nc subcarriers), and represent eacfi subchannel via an 
average gain. By doing so, we can treat the subchannel gains as being 
independent of each other. 



Theorem 2. Suppose that F{x, r) : R" x E"- R is a 
function of li. ^ M" and r € M"'', and r is a random vector 
whose components are nonnegative. For every e > 0, if there 
exists an yi ^ R" such that 

inf {vl/(x,^.)-j?e}<0, (25) 

where 

^-Cx, ^?) ^f?E{cxp(e-iF(x,r))}, 

then Pr {F(x, r) > 0} < e. 

Proof: (Sketch) The proof of the above theorem is given 
in IIT3I in details. To help the readers to better understand the 
idea, we give an overview of the proof here. 

It is shown in |fT3) (see section 2.2 therein) that the 
probability Pr{i^(x, r) > 0} can be bounded as follows: 

Pr{i^(x,r) > 0} < E{i/;(£»-iF(x,r))} . 

Here, > is arbitrary, and :R^-Ris a nonnegative, 
nondecreasing, convex function satisfying V'(O) = 1 
V'(-z) > V'(O) for any z > 0. One such ijj is the exponential 
function ip{z) = exp(z). If there exists a ^ > such that 

E{exp(^"^F(x,r))} < e, 

then Pr{i^(x, r) > 0} < e. By multiplying by ^ > on 
both sides, we obtain the following sufficient condition for 
the chance constraint Pr {-F'(x, r) > 0} < e to hold: 

*(x,^)-^e<0. (26) 

In fact, condition ( |26] | is equivalent to (l25l l. Thus, the latter 
provides a conservative approximation of the chance con- 
straint. □ 
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